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We develop a unified theoretical framework for the anisotropic Kondo model and the 
boundary sine-Gordon model. They are both boundary integrable quantum field theories 
with a quantum-group spin at the boundary which takes values, respectively, in standard 
or cyclic representations of the quantum group SU{2)q. 

This unification is powerful, and allows us to find new results for both models. For the 
anisotropic Kondo problem, we find exact expressions (in the presence of a magnetic field) 
for all the coefficients in the "Anderson- Yuval" perturbative expansion. Our expressions 
hold initially in the very anisotropic regime, but we show how to continue them beyond 
the Toulouse point all the way to the isotropic point using an analog of dimensional reg- 
ularization. The analytic structure is transparent, involving only simple poles which we 
determine exactly, together with their residues. For the boundary sine-Gordon model, 
which describes an impurity in a Luttinger liquid, we find the non-equilibrium conduc- 
tance for all values of the Luttinger coupling. This is an intricate computation because 
the voltage operator and the boundary scattering do not commute with each other. 
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1. Introduction 

One-dimensional quantum field tlieories witli gapless bullc excitations and boundary 
interactions display a wide range of interesting characteristics. They exhibit crossovers 
between Fermi liquid and non-Fermi liquid behavior, they can be sucessfuUy treated by a 
variety of powerful and interesting techniques, and they can be realized experimentally. 

The classic example of such a system consists of electrons interacting with dilute 
impurities in a metal, which can be described by the Kondo model. This system is actually 
three-dimensional; it can be described by a one-dimensional model because with dilute 
enough impurities, the interesting physics occurs in s waves around each impurity, and one 
can restrict attention to the radial coordinate. This model has a variety of experimental 
realizations, and has been the focus of much attention in the last 30 years (see and 
references within). Of recent interest is the problem of an impurity in a Luttinger liquid 
1^. A Luttinger liquid (interacting electrons in one dimension) may be realized in a one- 
dimensional wire or by the edge of a fractional quantum Hall device [|]. A fractional 
quantum Hall device is made by putting an electron gas trapped in two dimensions into 
a strong transverse magnetic field. When the Hall conductivity is locked to its plateau 
value, the current flows only along the edges of the device, and the system can be described 
effectively by a one-dimensional theory. Experiments have been done on the conductance 
through a point contact (which is the impurity in the theory) in one of these devices 0, 
and they agree well with theory 

The objects of our attention in this paper are one-dimensional models with interaction 
on the boundary only. We concentrate on two such models, the one-channel Kondo model 
and the massless boundary sine-Gordon model. The problem of an impurity in a Luttinger 
liquid can be mapped onto the latter. Moreover, when the bulk degrees of freedom are 
integrated out, both describe problems in dissipative quantum mechanics f^: a particle 
moving in a double well for Kondo (an infinite number of wells for boundary sine-Gordon) 
with a dissipative environment. 

In this work, we show that the Kondo model and the boundary sine-Gordon model 
can be treated in the same theoretical framework. Both can be reformulated as a free 
boson on the half-line interacting with a spin on the boundary, where the spin is in a 
representation of the "quantum-group" algebra SU{2)q. This algebra, as we will discuss 
below, is a one-parameter deformation of the ordinary SU{2) algebra. In the Kondo model 
the spin is in a standard spin-j representation, while for the boundary sine-Gordon model 
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the spin is in a "cyclic" representation, a quantum-group representation which has no 
analog in ordinary SU{2). We will find a simple relation between the partition functions 
of the two models. The relation is established through the use of the trace of the quantum 
monodromy operator, an object generating the conserved charges of the quantum KdV 
system Having this relation, we can relate quantities in one model to quantities in 
the other model. For example, the perturbative coefficients of the partition function of 
the spin— 1/2 Kondo model are expressed in terms of ordered integrals which are difficult 
to evaluate. This relation yields an expression for these coefficients in terms of known 
coefficients of the boundary sine-Gordon model . 



The starting point of the theoretical analysis is a one-dimensional quantum theory at 
a fixed point of the renormalization group. This means that the system has no mass scale, 
so there is no gap in the spectrum. Such a model can be described by a 1 -|- 1-dimensional 
massless quantum field theory. The issue of the boundary conditions in these models is not 
a nuisance but in fact can be of crucial importance. Basically, most of the physics which can 
happen in the bulk can also happen on the boundary alone. Studying boundary behavior 
is not only simpler mathematically, but it can also be easier to observe experimentally. 
There are very few experimental probes of one-dimensional quantum systems, and the 
ones mentioned above are both boundary effects. 

A boundary fixed point is a point where the boundary condition does not destroy 
the scale invariance of the bulk; the methods of boundary conformal field theory are 
applicable here. However, an interacting boundary condition as in both the above systems 
will introduce a scale to the problem, which we generically call Tb (in the Kondo problem 
this is often referred to as the Kondo temperature Tk)- Although by definition bulk 
effects in these models do not depend on this scale, boundary effects of a system at non- 
zero temperature can now depend on the dimensionless parameter T/Tb- Varying this 
parameter allows one to interpolate between different boundary fixed points. For example, 
in the Kondo problem at T/Tb oo, there is a boundary fixed point where the impurity 
decouples. As T/Tb 0, one approaches another boundary fixed point where the electrons 
bind to the impurity. (The properties of the low-temperature fixed point are far from 
obvious; it took years of effort to establish them.) For the boundary-sine Gordon model, 
the fixed points correspond to Dirichlet and Neumann boundary conditions on a boson; 
which is the high-temperature one and which is the low-temperature one depends on the 
the boundary coupling. 
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The field theories we discuss have the special property that they are integrable, as 
are many one-dimensional theories. As a result, we can do many calcinations exactly. In 
this paper, we mainly discuss the partition function and free energy. However, transport 
properties (which are experimentally measurable) can also be computed exactly . The 
methods we will describe enable one to study these systems for all values of the coupling 
— near and far from the fixed points. Other methods generally rely on perturbation 
theory around these fixed points. Another advantage, for example, is that in the Luttinger 
problem one can compute transport properties like the conductance even out of equilibrium. 
Standard field theory techniques are not applicable; the best one can do is use Kubo formula 
to calculate the linear response near V = 0. 

The models discussed here can be treated as a boson where the one-dimensional space 
is a half-line. The only interactions take place on the boundary. The bulk Hamiltonian 
takes the form 

Ho = -— da [U^ + {d^cPf] + — / dad^<p. (1.1) 

The first model we discuss is the boundary sine-Gordon model. The boundary hamiltonian 
is 

Hbsg = '2v cos (PiO). (1.2) 

In this case, the parameter V of ( pTTP plays the role of a physical voltage, while v is related 
to the boundary scale Tb in a manner to be discussed below. 

The second model, the one-channel anisotropic Kondo problem of spin j/2, can also 
be expressed in this form using the well-known technique of bosonization [jl2| . We ignore 
the charge sector of the Hamiltonian, which does not interact with the spin and decouples 
from the problem. The total Hamiltonian is then H = Hq + Hj, where the boundary 
interaction Hj= J2i=x y z ^i'^i^i- Here, Si is the impurity spin on the boundary, Ji are 
the fermion currents and li are the coupling constants. The problem is anisotropic when 
Iz ^ Ix = ly In the bosonized language, the boundary Hamiltonian becomes 

Hj = A (5+e^'^(°) + 5_e-^^(°)) - ^5,. (1.3) 

We have replaced the original parameters = ly and Iz with g and A; g parametrizes 
the anisotropy [g = 1 is the isotropic case and (7 = 1/2 is called the Toulouse limit) while 
A oc \Ix\- The precise relation is not universal so we will not need it here. The JzSz term 
has been absorbed in a redefinition of g. Traditionally in the Kondo problem, one takes 
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the matrices Si to act in the spin-j/2 representation of SU{2) (in (L3) and in aU what 
foUows we use the conventions that eigenvalues of Sz are integer, so they e.g. take values 
Sz = ±1 in the spin 1/2 representation). However, for the problem to be integrable, one 
must instead take them to act in the spin j /2 representation of the quantum group SU (2)^, 

[Sz, S±] = ±2S±, [S+, S-] = (1-4) 

where q = e^'^^ . In the isotropic case q = —1, this reduces to the usual SU{2) algebra. 
The distinction between SU{2) and SU{2)q is not important for j = 1 or j = 2 at ar- 
bitrary q, because the spin 1/2 representation remains the Pauli matrices and the spin-1 
representation is also the same up to a rescaling of Ss^ and S- . In the following, the Kondo 
model of spin-j/2 is defined to be the model with the g-deformed algebra, so it can be 
identified with the "physical" Kondo model only for j = 1, 2. In the Kondo models, dacj) is 
the 2;-component of the fermion current, so V corresponds to an external magnetic field. 

We consider the system in imaginary time compactified on a circle of length 1/T with 
T the temperature. Defining the partition function via the trace Zj = Tr 6*^^°+^^)/^, we 
introduce Zj = {j + l)Zj{X)/Zj{0) and Zbsg = 2^bsgW/2^bsg{0)- In the following, we 
often use the variable p defined as 

= 27rp. (1.5) 

Strictly speaking, the boson hamiltonians at non-zero temperature make sense only when p 
an integer. The consideration of real (physical) voltage or magnetic field requires analytic 
continuation. 

The dimension of the vertex operators e^^^ is g. For example the two-point function 
on the boundary is 

^g^^(o,r)g-^^(o,r')^ ^ ^ sin7rT(r - t') , (1.6) 

with K the frequency cutoff arising from the normal-ordering of the operators e^^^. We 
will denote the case l/2<(7<las the repulsive regime and < g < 1/2 the attractive 
regime; "attractive" and "repulsive" are the corresponding type of fermion interactions 
when one fermionizes this model into the Luttinger mo del 0. The Toulouse limit g = 1/2 



^ In the quantum wire problem (an impurity in a Luttinger liquid) where one starts with 
electrons on a full line, the entire domain < g < 1 corresponds to repulsive interactions between 
the physical electrons. There is a rescaling of the coupling when one maps the model on to the 
half- line M. 
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corresponds, of course, to free fermions. For g > 1 the vertex operators are irrelevant, and 
the model is best approached by using a "dual" picture 0. 

The paper is organized as follow. In section 2 the attractive regime is described. 
Results are obtained to all orders in perturbation theory using Jack symmetric functions for 
the boundary sine-Gordon model. Then, making use of the monodromy matrix, a relation 
is given between this latter model and the Kondo models. In the last part of this section, 
the thermodynamic Bethe ansatz is used to provide non-perturbative results in both cases. 
In section 3, the repulsive regime is explored. There the perturbative coefficients of the 
partition function diverge and a regularization is needed, which usually is provided by a 
high-frequency or short-distance cut-off. We show, using the explicit expressions discussed 
in section 2, that these divergences can also be controlled by analytic continuation from 
the repulsive regime, an analog of dimensional regularization. Coefficients for the free 
energy can be obtained in this fashion all the way to (7 = 1; at particular values of g 
there are poles, and we compute the residues exactly. These results are in agreement with 
computations using the Bethe ansatz in the repulsive regime, where the poles result in 
logarithmic terms in the free energy. In section 4, the relation between models is extended 
to non-zero V. This gives the perturbative coefficients in Kondo model as a function of 
magnetic field. Moreover, it yields the previously-unknown conductance for the boundary 
sine-Gordon model at all values of g. Some final remarks are collected in the conclusion. 



2. The attractive regime at zero voltage 

In this section we review earlier results for the anisotropic Kondo problem at zero 
magnetic field and the boundary sine-Gordon (BSG) model at zero voltage and with g < 
1/2. There are three useful and complementary approaches, all of which we will later 
extend to finite magnetic field (resp. finite voltage) and to (7 > 1/2. 

We first discuss how to expand the partition function in powers of the interaction 
strength. For the Kondo problem, this was ffist considered long ago in |1^,12|, where the co- 
efficients of this expansion were expressed as multiple integrals. These integrals are rather 
complicated, and until now had not been evaluated explicitly except in very special limits. 
The partition function in this form is equivalent to that of a one-dimensional gas of posi- 
tive and negative charges with logarithmic interactions (equivalently of a two-dimensional 
Coulomb gas on a circle). For the boundary sine-Gordon model, the perturbative expan- 
sion is formally very similar, but not identical. In that case, the multi-dimensional integrals 
can be explicitly evaluated, using recent results for symmetric polynomials |]1^ . 
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The second approach uses integrabihty, albeit in a rather abstract way. We define 
the trace of the "quantum monodromy operator" 0, whose expectation value gives the 
Kondo or BSG partition function, depending on which representation is chosen |T^,^|. 



Using some properties of this operator, we are able to relate the Kondo partition function 
to the BSG one. As a result we are able to evaluate the integrals in the Kondo expansion 



explicitly, using the already-evaluated BSG ones fTj,^ . 

The third approach uses integrabihty in a more standard way. We describe the model 
in terms of interacting quasiparticles and their scattering matrices. The thermodynamic 
Bethe ansatz (TBA) can then be used to derive the free energy and related quantities for the 
Kondo model and for the BSG model |jT^. The direct relation between the partition 



functions can be rederived, at least for values of the coupling g = 1/t, t integer. The 
TBA approach has the disadvantage that at non-zero temperature the integral equations 
derived are not continuous in g (although the final results of course are). However, it has 
the advantage that it allows transport properties like the current and conductance [[7i|TT| 



and the zero-temperature noise [jT^ to be computed for the BSG model. (For Kondo, only 
the zero-temperature magnetoresistance has been computed Some simple relations 

have also been derived relating transport properties to equilibrium properties |T^; we 
generalize these in sect. 4. 

2.1. Perturbative approach 

The partition functions Zj and Zbsg can be expanded in powers of A and v respec- 
tively. The term of order A^" or f involves a correlation function of n vertex operators 
e*"^ and n vertex operators e"*"^, all living on the boundary. These multi-point functions, 
evaluated in the free-boson theory and by Wick's theorem, are reduced to a product of 
two-point functions like ( |1.6|) (see e.g. for a review). The problem then becomes for- 



mally equivalent to a two-dimensional Coulomb gas with positive and negative charges 
restricted to live on a one-dimensional circle. To calculate the partition function, we must 
integrate over the locations of the charges. The integrand is then the scaled correlator 



2^2n(K}, {W-}) 



^ /J(t}iui/27TT) z<^(u„/27rT) -i</)(u;/27rT) 

2nT) ^ 



n.<,4sin(^)sin( 



/ I \ 2q 



n,,2sin(^ 



f2.1) 
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The difference between the Kondo model and the boundary sine-Gordon (BSG) model lies 
in the limits of integration. In the Kondo model, each of the vertex operators comes with 
a spin operator; the thermal average of monomials of vertex operators are computed as 
in ( |2.1|) , while one has to take the trace of the corresponding monomial of spin operators 
in the representation of interest. This puts various contraints on the order of the charges. 
For example, S'^ = S'^ = when the spin is 1/2, so only terms of the form Sj^S-Sj^S- . . . 
survive in the perturbative expansion, and consequently charges alternate in sign on the 
circle. Thus in terms of the renormalized parameter x, defined as 



the spin- 1/2 Kondo partition function is [T^,|TB 



Z^ix)=2 + J2^^'"Q2n, (2.2) 



where 

r"27r fUi 



Q2n{p) = 2 / dui / du[ / dU2... du'^l2n{{Ui},{u'i}). (2.3) 

Jo Jo Jo Jo 
The effect of the charge ordering is seen in the limits of integration. Higher-spin partition 
functions have the same integrand, but with the appropriate resctrictions on charge order- 
ing. For the boundary sine-Gordon model there is no boundary spin, so one has unordered 
integrals: 

oo 

ZBSG{x)=l + Y,xl^SGl2n, (2.4) 



where 



and 



n=l 



XBSG = ^{ — 



I2n = -r-rr2 dui... du'^J2n{{ui},{u'i}). (2.5) 
l^'J Jo Jo 

The lack of ordering makes no difference for n — 1, so I2 = Q2j but the others are different 
I. 



^ The l2n were denoted Z2n in [10,14|; we change notation here to avoid confusion with the 
higher-spin partition functions. 
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The unordered integrals In can be computed exactly in terms of an n-dimensional 
series |T^: 

T 1 v^TT / r[m, +^(7i-z + l)] y 

where the sum is over all sets (Young tableaux) m = (mi, m2, . . . , m^), with integers 
rrii obeying mi > m2 > ... m^ > 0. For n = 1 the series can be summed, giving 
I2 = r(l — 2g)/[r{l — g)]"^. Although this series looks quite imposing, it can be generated 
by a simple recursion relation. We introduce the truncated sum /2ri(A), which is defined 
as the sum over m with the condition that all m^ < A. Then it is not difficult to show 
that 

These relations allow a precise determination of the partition functions up to large orders 
in the perturbation expansion. 

The boundary sine-Gordon model can be placed in the same framework as the 
anisotropic Kondo models. It is enough to discuss the simplest case when g is a root 
of unity {g rational), since by continuity, the results we derive will hold for any q of unit 
modulus. Suppose therefore q'^ = ±1 for some integer k, and consider a cyclic representa- 
tion of the quantum group SU{2)q (see [0). These representations, which are labeled by 
an arbitrary complex parameter 6, have no highest- or lowest-weight state; the states are 
eigenstates of Sj^ or S~ to the tth power. They have dimension k, with a basis of states 
|m > such that 

q{S-m)/2 _ ^-{S-m)/2 

S+lm > = — |m -|- 1 > 

q-q 

qiS+7n)/2 _ ^-{S+m)/2 /2.8) 

S-lm >= — |m — 1 > 

q-q 

Sz\m >=m |m >, 

where states |m > and |m mod k > are identified, and the fundamental set is chosen to 
be 0, 1, . . . , (/c — 1). To obtain the BSG model, we set q^ = C and let C » 1 and real (so 
5 is imaginary). Thus 

=F"^/2 

S±\m>^C- -|m±l>, (2.9) 

q - 

so in this limit the commutator of S± can be neglected and the traces of all monomials 
become identical: 

/ 1/2 \2n 

TrM = k ^ _^ , (2.10) 
\q-q V 



where Ai is the product of n operators S+ and n operators 5"- in any order. Thus when 
evaluating Zg for C large, all the possible orderings of S+ and S- within the trace have 
the same weight, so 

Zs{x) ^ kZssG (c ^ ^ , x] , / = C » 1. (2.11) 

This observation allows us, for example, to find the boundary S matrix of the BSG model 
pod , as we detail in the Appendix. It will also enable us to derive many properties of the 
partition function Zbsg in the subsequent sections. 

2.2. Quantum monodromy and fusion 

We introduce the quantum monodromy operators associated with these models M 



Lj{x) = Uj <( e'^^^^Pexp 



(2.12) 

where Ilj indicates that the matrices Sj are in the spin-j/2 representation, V indicates path 
ordering, and the exponentials are normal-ordered. In this formula, P is the momentum 
operator appearing in the mode expansion of the left- moving field (pL 



^L(r) = Q + 27vTPt + z 2^ -f e 



O-n -2iiTnTT 



n 



Here the normalizations are [(5,P] = ig, [a„, a^,] = ^Sn+m- Observe that Lj is an 
operator acting both on the spin degrees of freedom and on the "free-boson" degrees 
of freedom. By expanding Lj in powers of x and noting that with Neumann boundary 
conditions 2(^^(0, r) = ^(0,r), one finds that the partition functions are equal to the 
eigenvalues of the quantum transfer matrices acting on momentum eigenstates P\p) = p\p), 

Zj{x,p) =< p\tr e'''^^^Lj{x)\p >, (2.13) 

where the trace is computed over the spin degrees of freedom. At zero voltage, p = 0; we 
define Zj{x) = Zj{x,p = 0). 

As observed in [0, the Lj satisfy the Yang-Baxter equation. Using the fusion of 
quantum transfer matrices, one can prove the identities 

Zj{q^/^x)Zj{q-^/'^x) = 1 + Zj-i{x)Zj+i{x) 

Zi(g(^+i)/2x)Z,(x) = Zj+iiq'/^x) + Z,_i(q-V2^) (2.14) 
Z^{q('+'y'x)Zs{x) = Zs+i{q'/'x) + Zs-i{q-'/'x). 
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The first was discussed in ; the second can by proven by using the first and by induction. 
The last foUows using the same technique as in | , together with the fusion rules for cyclic 
and standard representations ||2^ 



Using these relations together with ( |2.11| ), we can express the boundary sine-Gordon 
model partition function in terms of the Kondo partition function. We have, from ( |2.14|) , 

Zx{Cq^l'^x)ZBSG ( — — tx] = Zbsg (— — tx] + Zbsg (— tx] 

\q-q 1 / \q-q ^ J \q - q ^ J 



from which it follows that [|n|,|T5 



\(a - a-^)T] - ZBSG{qx) + ZBSG{q 'x) 
Z,[{q q )x]- ^--^^ . (2.15) 

Inserting the perturbative expansions into ( |2.15| ) gives the Q2n in terms of the already- 
known thus completing the derivation of the perturbative partition function for g < 
1/2. 

2.3. The thermodynamic Bethe ansatz 

The fusion relations discussed in the previous subsection are one of the many conse- 
quences of integrability ||2^. The standard way of approaching the problem is to use the 



Bethe ansatz. Here, one derives integral equations which determine a set of functions ej{9), 
where ^ is a rapidity (the logarithm of the energy of an individual particle). The ej{9) can 
be thought of as the energy of an interacting quasiparticle, in the sense that the energy 
of the entire system shifts by Tej{9) when a particle of rapidity 9 is added to the system. 
Moreover, the distribution function is given by 1/(1 + exp(ej)). Many physical quantities 
can be expressed in terms of these functions. Since this approach has been discussed in 
detail in many places, we start with the integral equations and discuss their consequences. 
For technical reasons, we consider the case g = 1/t, where t is an integer. The integral 
equations for both Kondo and BSG are [l],0,|lf 



where the "incidence matrix" Njk is defined by the diagram 



12 k 

o — o o- - -o — O t 
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+ 



where Njk = 1 if the nodes j and k are connected, and Nj^ = if not. The solution of 
these integral equations is fixed uniquely by demanding the asymptotic form 

ej ^ 2 sin e± ^ as ^ ^ oo. (2.17) 

z[t — 1) 

This asymptotic form is just the energy of the individual particle over T; the interactions 
become negligible in the large-energy limit (equivalent to sending A and v to zero in the 
Hamiltonian) . 



The free energies in this regime for Kondo (with spin less than t/2) and BSG 
can be written in the form 

^ _ sin[j7r/2(t - I)] ^ f dO t-1 

^^■=^^cos[7r/2(t-l)] 2^cos\y[{t-l){e-\uTs/T)] ^^(1 + ^^)' ^^.IS) 

for J = 1 ... t — 2 together with Ft-i = 2Fbsg and 

Fbsg = :: TTf^-nT-^/ o vjT;— n?^— nFT^ + ' ^2.19) 

2 cos[7r/2(t — IjJ J 27r cosh[(t — 1}(6^ — InTB/rjJ 

where et-i = €+ = €-. 

It has been shown that the equations ( |2.16| ) require that ej{9) = €j(9 + z27rt/ (t — 1)) 



||24|| . This means that the integrals in ( |2.18| ) and (|2.19| ) can be expanded as a power series 
in (Tb/T)^*-*"^-*/*, so we see that the bare couplings A and v and the renormalized coupling 

1-9 



X are proportional to . In fact, for BSG, the exact constant was determined in |jTT[ 

and is 

for any value of g, not just g = 1/t. At fixed T^, the Kondo bare coupling A is related to 
the bare BSG coupling via a constant to be determined at the end of this section: X = ^v. 
This constant ^ is independent of the impurity spin considered, as observed in 0. With 
this relation of x and Tg, we see that the second term in (|2.18| ) or ( |2.19| ) is an analytic 



power series in x, like the perturbative partition functions in sec. 2.1. 

We must take care in relating these non-perturbative free energies to the perturbative 
partition functions discussed before. The TBA deals with excitations over the vacuum; 
by convention, the ground state (no particles) is assigned a vanishing energy and entropy. 
Therefore, one expects Fj and Fbsg to equal to the perturbative partition functions defined 
previously up to a constant shift (which on dimensional grounds must be proportional to 
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Tb) and a term proportional to T. Neither of these change e.g. the specific heat. This 
ambiguity is fixed by studying the behavior at Tb = and by studying the analyticity 
properties. From the TBA equations, it is simple to derive from ( |2.16| ) that as ^ ^ — oo, 
the functions ej{9) go to a constant, which is 

e^^(-°°) = (i + 1)2-1, e^±(-°°) =t-l. (2.21) 



Plugging this into (|2.18| ) gives 



FjiTB = 0) = -Tln(l +j), Fbsg{Tb = 0) = -|lnt. 

This fixes the piece proportional to T. The piece proportional to Tb is fixed by noticing that 
because the perturbative partition function is analytic in x oc ^, the term proportional 
to Tb cannot appear here. Thus the relation between the perturbative partition functions 
and the TBA free energies is 

^ _ sin[j7r/2(t - 1)1 
Fj = -TlnZj + Tb- ^ ^ 



cosK/2(i-l)] 

Tb T ^ ^ 

FbSG = -TIuZbsG + 75 r^TTTT TTT " TT^^^ 

2 cos[7r/2(t — 1)J 2 

The role of the shift is simple: it precisely cancels the large Tb/T behavior of the partition 
functions. As is easily seen by substituting the asymptotic form ( p.l7|) into ( p.l8|) and 
( P.19|) , the free energies Fj/T and Fbsg/T determined by the TBA go to zero as T /Tb 0. 
Meanwhile it was shown in [l^] that — In Zbsg oc Tb/T in this limit, and the relations 
( ^.15| ) and ( ^.141 ) indicate that — In Zj grows as well. In fact from ( |2.22D it follows that 



Zj ^exp 



Tb sin[j7r/2(t - 1) 
Tcos[7r/2(t- 1)] 



^ I f Tb 

Zbsg exp 



^/t \2Tcos[n/2{t-l) 

Thus we see that although the perturbative partition functions grow exponentially for 
large x, the series expressions are still convergent (they actually have an infinite radius of 
convergence) . 

Since we have related the TBA results to the perturbative ones, we can combine the the 
fusion relations ( p.l4| ) with the TBA equations ( |2.16| ) to give much more information. For 
example, the relation ( |2.18| ) gives the Kondo partition functions Zj only for j = 1 ... t — 2, 
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but the remainder can be generated from (|2.14 ) . The relation ( 2.22| ) allows us to write the 
perturbative partition functions in terms of the ej very simply. Denoting convolution by 



—A{a-9)B{9), 
-oo 27r 



one has 

lnZj{x) =st-i *ln(l + e'0(") j = l...t-2 
InZBSciCx) =-^lnt + st-i * ln(l + e'^-'){a) (2.23) 
lnZt-i{x) = 2st-i * ln(l + e^*-i)(a), 
where Sa = a/ cosh(a^) and a = In Tg/T. These relations also give immediately 

Zt-i{x)=tZlsGiCx)- (2.24) 

In the following, we often switch from the x variable to a, keeping these relations in mind. 
Introduce as usual 

Yj{a) = e'^^''\ (2.25) 

The Yj are analytic functions of oc exp[2ta/(t — 1)] |2^]. From the TBA equation for 
ei one finds then Yi{a) = Z^ix) and from the TBA relation for et_i, Yi_x{a) = Zt_2(x). 
The TBA relations for the other nodes then imply that 

Y,{a) = Z,+i(x)Z,_i(x), i = 1, . . . ,t - 2. (2.26) 

This simple relation between the TBA and the perturbative partition function gives exact 
perturbative expressions for all the ej. This is consistent with the relation ( |2.24| ), since the 



original integral equations ( |2.16| ) along with ( p.23|) give Yt-i = tZtsZ^gQ 



2 

We can in fact rederive the results of sect. 2.2 at q = e^'^l^ by converting the TBA 



equations (|2.16| ) into functional equations in the complex a plane. Using the identity 



Sa(^ + ^7^) + «a(^ - ^7^) = 27r5(0). (2.27) 
la 2a 

with {^^) gives 

Zj{q^/^x)Zj{q-^/^x) = 1 + Yj{a). (2.28) 



where j = 1 . . .t — 2. Plugging the relation (|2.26| ) into ( |2.28|) recovers the fusion relation 

(EH) 

Zj{q^/^x)Zj{q-^/^x) = 1 + Zj_i(x)Z,+i(x) (2.29) 
13 



from 0. Once written in terms of the variable q, (|2.29| ) holds for q generic, as discussed 
in sect. 2.2. Observe that such a direct proof of ( |2.29D establishes conversely the fact that 
the relation between Tg/T and x is independent of spin. 

When q = e^^/*, t an integer, these fusion relations close i. We have from (|2.29|) 



Zt-iiq^/^x)Zt-iiq-'/^x) = 1 + Zt-2{x)Ztix) 



while using (|2.27| ) in ( |2.23|) gives 



ZBSG{^q'^^x)ZBSG{^q-'^^x) = -{1 + Yt-i{a)). 



(2.30) 



Thanks to the identification Yt-i = Zt-2 these two relations are compatible with (|2.24| ) if 
and only if 

= Zt_2(x) + 2. (2.31) 



The latter relation follows from the quantum group representation. Indeed, when q is 
a t^^ root of unity, the rep 
schematically as in figure 1. 



a t^^ root of unity, the representation of spin t is reducible because 5"^ = 0, and looks 



m=t-1 



m=t-2 



m=t-3 



m=0 



f 



Fig. 1: A schematic representation of a cyclic representation of SU{2)q 
when g is a t^^ root of unity. Up and down arrows represent the action of 
raising and lowering generators. 



^ A different closure happens in the minimal models of conformal field theory Q, where the 
end nodes t — 2, + and — are removed from the incidence diagram. 
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We see that the states with values = ±t do not contribute to the trace of any 
monomial in Sj^S- of non vanishing order. Hence in the perturbative expansion, all terms 
for spin t — 2 and spin t are equal, except the term of order zero that simply counts the 
number of states. This term differs by two in the two representations, and ( |2.31| ) follows. 

We can also rederive the relation ( |2.15| ) between the Kondo and BSG partition func- 
tions without using cyclic representations by using the fusion relation from ( |2.14| ): 

Zi{ix)Zt-i{x) - Zt-2{q-^'''x)+Zt{q^'^x). 



Using (|2.31D , we rewrite the right hand side as 

1 + Zt_2(?^/'x) + 1 + Zt-2{q-^'^x) 



which using Yt-i = Zt-2 and (|2.3UD this is in turn 



tZBSGi^x)[ZBSGi^qx) + ZBsoi^q ^^)] 
from which, using (|]2|), it follows that 

[fa + (2.32) 

Zbsg{x) 

together with the fact that ^ = . This value for ^ is obtained simply by matching 

the first order in the perturbation theory, since we know that I2 = Qi- 



3. The repulsive regime at zero voltage 

3.1. Poles and log terms in the perturbative expansion 

The previous section concerned the attractive regime g < 1/2. This, for example, 
is the regime of greatest interest in dissipative quantum mechanics, where the particle 
exhibits oscillatory behavior. The filling fractions z/ = 1/(2?! + 1) where the edge modes in 
the fractional quantum Hall effect are described by the BSG model also lie in this regime. 
However, the original isotropic Kondo model is at (7 = 1, and we will see that there is a 
great deal of interesting behavior in the repulsive regime 1/2 < g < 1. We will address 
this regime largely by exploiting some simple analyticity properties. In particular, we will 
show how to obtain an analytic expression for the coefficients all the way to g = 1. 
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We will show in this section that if we define the expansion of the free energy of the 
spin- 1/2 Kondo model as 



-T In Z^=TJ2 f2nX^" = -T ln(2 + J] Qsr 

n=0 n=l 

the coefficient f2n has a simple pole aX g = 1 — l/(2n), with residue r2n = — l/(27rn^). The 
half-integer-spin Kondo and BSG free energy expansions have poles in the same places. 
Moreover, when this divergence is regulated properly, we find a term —2nr2nTB ^^^{Tb/T) 
in the free energy for g = 1 — 1/ (2n). This term yields, for example, a term linear in Tb/T 
in the specific heat, indicating that these values of g (which include the Toulouse limit 
g = 1/2) are pathological in some respects @. 



The first thing to notice is that the integrals ( |2.3|) and (|2.5| ) , which define the pertur- 
bative coefficients l2n and Q2m diverge at short distances when g > 1/2; correspondingly, 
the series expansion ( p.6|) diverges for g > 1/2. There are a variety of ways to regulate the 
integrals. In a numerical approach, this would be done using a cut-off. However, the most 
natural approach here is analytic continuation. This approach, which is very analogous to 
dimensional regularization [^, means we define the regularized integrals as the analytic 
continuation of their values for g < 1/2. This continuation is illustrated by examing the 
first coefficient, /2 = —/2/2, which we saw in sect. 2.1 is given by /2 = r{l — 2g)/[r{l—g)]'^. 
At (7 = 1/2, /2 has a simple pole. There are no branch points anywhere, and since it is 
finite for all other g < 1, the analytic continuation is perfectly well-defined. Implicit in the 
following is the assumption that the regularization done for the Bethe ansatz (the cutoff 
of the Fermi sea for Kondo 0) gives the same results as this analytic continuation from 
g < 1/2. This assumption is certainly physically obvious, since by defining renormalized 
parameters one removes all cutoff dependence from the Bethe ansatz. Moreover, in the 
BSG model one starts directly from the regulated theory with no cutoff dependence [pO 



We now show that the large-Tg/T behavior of the partition function requires that 
f2n^^^ and f2n have simple poles aX g = 1 — l/{2n) for all n. As discussed in sect. 2.3, in 
this limit —TIuZbsg behaves like 

1 



-Tb- 



2 cos 



2(1-3) 



^ However, notice that now a term T'^"^^^^^ appears in the specific-heat expansion at all g; the 
log term is required to make this true at = 1 — 1 /2n. Thus in this sense the log terms make are 
not pathological but instead make the values g =1 — l/(2n) more like other values of g. 
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while —TliiZi goes as 



-Tb tan 



,2(l-£ 

Our analyticity assumption implies that these hold for (7 > 1/2 as well. Notice that these 
expressions have a simple pole as g 1 — l/(2n). In the TBA free energy this term is 
subtracted off, as seen in ( p.22|) . Because the TBA is finite, this divergence therefore is 
matched by one in the perturbative expansion. Since x oc {Tb/T)^~3 , the terms f2n^'^^x'^^ 
and f2nx'^^ are proportional to Tb when g = 1 — l/{2n). Therefore, there must be a simple 
pole in f2nX^^ at (7 = 1 — l/(2n), with residue = Tb / {'2.Ti'n?T). The pole in f2^^^^x^^ 
has residue (— l)""'"^r2n/2. By the same argument, the free energy coefficients in the spin- 
i/2 Kondo model when j is odd each have a single pole at (7 = 1 — l/(2n) with residues 
r2n- The free energy coefficients for the integer-spin Kondo model have no pole at these 
values. 



We can see these poles explicitly by studying the series expansion (|2.6| ) for the l2n in 
the boundary sine-Gordon model. Initially the series looks useless, because it diverges for 
g > 1/2, and we only know how to resum it for I2. However, a first interesting observation 
is that the series expressions for the f2n^'^^ converge even where those for the individual 
l2n do not. This is because some of the divergences in the Coulomb integrals are cancelled 
when taking the connected part. More precisely, we define the truncated series /2n(A) as 
the expression (|2.6|) with all < A and 



/f^^^(A) = 



(/2(A)) 



2 



- /4(A) 



Then, using the previously obtained recurrence relation (p. 71) , one finds that 

/i"'°'(A) - tr^HA - 1) ^ (^^0±l"^ r(g)-A^.- (3.1) 

for A large. This expression converges as A ^ 00 for g < 3/4. Moreover, the pole 
at g = 3/4 is clearly identified, and its residue can easily be computed, because the 
divergence is proportional to that of the zeta function. One confirms the earlier result that 
near g = 3/4 



ABSG) A _ 1 

167r(^7 - 3/4) ^ T 

where we used the relation ( p.20|) to relate x and Tb- 
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4(1-3) 



Since (|3.1|) tells us explicitly how the series diverges, the continuation around the pole 
can be constructed by adding and subtracting a zeta function. More precisely, we define 
the continuation to be 

ABSG), X ABSG)reg, x , 25f(l - 25f) + 1 

where 

.(B5G)reg. .X ^ ABSG) , f.. _ 2g(l - 2^) + 1 ^ 4^.4 
/4 (AJ_/4 (AJ 2(l-2^7)r(^7)4^^'' 



The "regular" part f^^^'^^'^^^ jg f^^^'^^ with the diverging part of the sum subtracted off. 
This series converges when A ^ oo for g < 1. One can extend this result past g = 1 in the 
same manner. 

It should be possible to find all the f2n for g < 1 in this manner. One first writes the 
higher f2nS in terms of the using the relation 



ABSG) _ ^ (-l)'(-)-^(/(m)-l)! 



(3.3) 



where m = {mi, m2, . . . TOi(m)} is a partition of an integer so that X] "^i = -^2(n) = 
-f2mi-^2m2 ■ ■ ■ ^ud Xj is the multiplicity of the integer j in m. We have checked that 
/f ^^)(A) - /f ^^^(A - 1) ^ C,{g)A^3-^ and /f ^^^(A) - /f ^^^(A - 1) ^ Csig)A'^-' 
when A is large, with CQ{g), Cs{g) known expressions. Thus poles in appear at 

g = 1 — l/2n for n = 1, 2, 3, 4 with the appropriate residue. We can then apply the same 
zeta-function method and regularize the sums to go all the way to g = 1. 

We have checked that the numerical values agree very well with the Bethe ansatz 
results. This takes some effort because finding the numbers from the Bethe ansatz requires 
that we numerically solve the integral equations, and then numerically fit the results to a 
power series. Moreover, at (7 = 1, we have Zbsg{x) = 1, We plot the results for fi 
and f^^^'^^ in figs. 2 and 3. 

We clearly see the pole in the data, and that f^^^^ and f^^^^^ do indeed go to zero 
as (7 1. Another interesting consequence is that it allows a very simple approximation 
formula for the f2n or f2n^'^^ for g near 1. We approximate the function by its pole 
plus a constant piece. For example, if we define the constant piece by requiring that 
f2n^^\9 = 1) = 0, we have 



2v^n2r(n)r(l - l/2n)2" 
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l-g-{l/2n) 



CD 
03 



^1 




-0.1 



-0.4 



Fig. 2: The free-energy coefficient f^^^^^ as a function of g. The pole is at 
9 = 3/4. 



C3 



0.15 



0.05 




-0.05 



-0.15 



-0.2 



Fig. 3: The free-energy coefficient f^^^^^ as a function of g. The pole is at 
9 = 5/6. 



The presence of these poles has interesting physical consequences. They indicate that 
the free energy at = 1 — 1/ (2n) cannot be expanded as a power series, but has an 
additional logarithmic term. The TBA free energy does not have a divergence, because 
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the pole is subtracted off, as in (|2.22| ). However, there is a leftover piece: 

2n(l-g)-l- 



lirrir. 



»l-l/(2n) 



Tb tan 



7rg 



2(1-^7) 



+ 



Tr 



2n 



g-l + l/{2n 
-2nTr2n In 



Tb 
T 



Tb 
T 



+ . 



Thus the free energy contains a logarithmic correction, arising from proper regularization 
of the divergence. Such terms are not unusual; for example they occur in the bulk free 
energy of the 2D Ising model and its (even) multicritical generalizations P^. As we will 
see in the next subsection, the existence of the log term also follows from the detailed TBA 
analysis. The fact that the free energy defined by analytic continuation has a simple pole 
at does not mean that the physical free energy - obtained with a cut-off regularization 
- diverges, but rather indicates that it has a logarithmic dependence on the cut-off at that 



point 



These results can be checked in the Toulouse limit g = 1/2, where the model is 
equivalent to a free fermion in a boundary magnetic field. As discussed in P] for example, 
the free energy for spin- 1/2 Kondo is 



-2T 



d9 



1 



27T cosh{6 - IuTb /T) 



ln(l + e" 



T 



B 



2Tlnr ( -j; 1 -2Tlnr 



■B 



2nT 



TT 



l-21n2-ln 



T 



B 



2kT 



Tin 2. 



Using the gamma function identity V{2a) = V{a)V{a + 1/2)2^" ^ I and the expansion 

lnr(a + 1/2) = lnr(l/2) + ^ ^^"-'^(1/2)^, 



n=l 



where ip^'^^x) is the mth derivative of the digamma function '(/^(x) = T'{x) /r(x), one finds 



-Tln2 + — (1 

TT 



ln^)+2Ty7/.("-iVl/2)-^ ^ . (3.4) 



Thus we see explicitly the log term oX g = 1/2, with coefficient = l/27r as derived above. 
Moreover, we see that all the /2n for n > 1 are finite at (7 = 1/2. This means, for example, 
that there is a double pole in Q4 at (7 = 1/2 in order for = —Q4/2 + Q2/8 to remain 
finite. In fact, this means that there is an nth order pole in Q2n aX g = 1/2 and that they 
are analytic in the neigborhood (and in fact all the way to (7 = 3/4). Moreover, aX g = 1/2, 
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Fbsg = and we have checked numericaUy that the values for from (|3.4|) are 

obtained by taking the hmit of the series expression as (7 ^ 1/2 ( p.6| ). 

A final comment is in order. At the isotropic point g = 1, the BSG model is trivial 
with these boundary conditions, so Zbsg{x) = 1 and l2n = 0. Notice that this follows 
easily from the relation ( ^.15| ). However, the Kondo problem is not trivial at g = 1 (this is 
the value of most physical interest), but the power series expression (proportional to T^"^) 
obviously requires modification. The fact that the exponent is vanishing is an obvious hint 
that there are log terms at every order, and indeed this is seen in the TEA solution |jl|,|l. 
Notice that the shift between the TEA and the power series has an essential singularity as 
(7 ^ 1, so if subtracted appropriately from the power series as g 1, the result may be 
finite and give the series with log terms at g = 1. We have not yet succeeded in carrying 
out this analysis. 

3.2. The TBA in the repulsive regime 

The TEA equations for the Kondo problem in the repulsive regime were derived in 
0. For technical simplicity, we consider only g = l — l/s,s>2an. integer. The equations 
are very similar to those in the attractive regime: 

where the incidence matrix Njk is as in sec. 2.3 with t replaced by s. The Kondo free 
energy is 

^^--^/f cosh(.-lnT,/r) '°('+-'") 

r M 1 ^^-^^ 

F,_i = -2T / — _ln(l + e-^-0- 

' J 27rcosh(^- In Ts/T) ^ ^ 

where e+ = e- = es-i. Even though the ESG problem is integrable in this regime, applying 
the TEA is difficult technically since now both the bulk and the boundary scattering 
matrices are not diagonal. We will use analytic continuation again to provide the ESG 
free energy. 

Defining this time 

Yj{a) = e-^^("), 
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we see right away that for j = 2 . . . s — 2, 

Yj =e-^^-+i/'^e-^-i/'^ (3.7) 

analogous to (|2.26[ ). Arguments identical to those in the attractive case require that 
y(a + zs7r) = Y{a) [^. Thus Y can be expressed an analytic power series in oc e^"/'^ = 
(Tb/T)2*^^~^) as before. Therefore, ( |3.7| ) indicates that for j even, Fj +e" is a power series 
in x'^. When j is odd and s is odd, Fj is a power series in as well, but for s even and j 
odd, any other term is allowed as well. In fact there is a log term as discussed in the last 
subsection. 

We can find the log terms at g = 1 — l/(2n) (i.e. s even) directly from the TBA, by 
plugging the power series expansion for Y into (|3.6D . For example, for s = 4 ((7 = 3/4), 
Yi{a) = 3 + ae"/2 + 6e°. Then, we see that 



F,/T + \n2-hx^ = - I f , \ ^ [ln(l + n(^))-ln4- V/^ 

J 2tx cosh[6 — mTB /T) L 4 

For Tb/T small, this is approximately 

de 1 ^^/4 + ae^/2 + (6-aV8)e^ 



27rcosh(^- In Tb/T) V 4 + 

c/6' 1 (6 - aV8)e^ 

2^ cosh(^- In Tb/T) 4 + 

4Tr „ /-"^ , w3 1 



(6-—) / du 



6-aV8TB, Tb 

'' \ — In 

47r T T 



One can in fact verify using functional relations analogous to those above and in [|T^ that 
h — a^/S = —2, so the coefficient is indeed — 4r4 = —Tb/{2'kT) as shown in the previous 
subsection. 

We can now derive the analogs of the fusion relations (|2.14| ). We define 

Zj(x) = e-^^^")/"^ jodd 
Zj{x) = e-^^(")/'^e-^" j even. 
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Using (|2.27|) , we have 



Z,((-Q)V\)Z,((-g)-i/2x) = 1 + 1^- = l + Z,_i(x)Z,+i(x) (3.8) 



analogous to (|2.28|) and ( p.29| ) in the attractive regime. The crucial difference is that q has 



been replaced by —q~^- The Zj (x) are analytic functions of for j odd, but for j even they 
include the log term, so implicit in this equation is the prescription —in < Imlny < in. 
Because the analytic continuation of Zj should still satisfy the fusion relation ( |2.29D , not 
all of the Zj can be the analytic continuation of the Zj from the attractive regime to the 
repulsive regime. However, notice that if we make the identification 



Tb sin jn{s - l)/2 
T cos7r(s-l)/2 



In z, (x) = In z, (x) + ^ ~ j odd 



In Z, (x) = In Z, {^x) + ^ 3(J- l)/2 ' 

then the Zj satisfy the fusion relations ( p.29|) . The shift as before cancels the pole in In Zj 
for j odd. Since Zj with j even is a power series in x^, the effect of the argument ix is 
to flip the sign of every other term. This is merely a matter of convention. With our 
choice q = e™^, ? = 1 for the classical limit (7 — while g = — 1 at the SU{2) point 
(7 = 1. Representations of SU(2)q and SU{2)_q-i are identical for j odd but they differ 
by a factor of i in the matrix elements of S± for j even. The coupling renormalization for 
j even would disappear if we chose to change the quantum group conventions. 

We can now flnd Zbsg by analytically continuing the functional relation (|2.15|) . Since 



this relation involves only q and the functions are series in x^, we can replace q by —q~^- 
With this replacement, all functional relations derived in in sect. 2.3 apply to the repulsive 
regime with t replaced by s. In particular, we showed that 2Fbsg{Cx) = Ft-i{x) implies 
( ^.15| ). Since the relation ( |2.15| ) for (7 7^ 1 determines all of the BSG coefficients 
uniquely in terms of the spin-1/2 ones Q2m given a Zi, it determines Zbsg uniquely. 
Therefore, we can reverse the argument in sect. 2.3, and say that given ( |2.15| ), we must 
have 

2FBSGi^x) = F^.iix) (3.9) 
where ^ = i/{q — q~^) as before. A similar result (missing the crucial factor of 2 and 



without specifying the ^) was conjectured in [^^. We emphasize that this relation is only 



true for s integer. The relation ( |2.15| ) of course is true for any value of g. However, the 
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exact result at integer s is very useful, allowing us for example to find the conductance in 
the BSG model exactly at these values, without any analytic continuation. 

We have checked the result (|3.9| ) numerically at s = 3 (again by comparing analytic 
continuation results to TBA ones) and find good agreement. The relation ( pj.20| ) still holds 
in the repulsive regime (the derivation of [|Tl|] holds for all g) , and we confirm also the value 



4. Non-zero voltage 

We now extend the results of sect. 2 to allow for non-zero voltage in the BSG problem 
and non-zero magnetic field in the Kondo model. For the Kondo problem, the TBA 
analysis is easily extended to non-zero magnetic field . The analysis is straightforward 
because the magnetic field couples to a conserved charge, the ^-component of the spin. 
Since the charge commutes with the Hamiltonian, the same diagonalization applies even 
with a magnetic field. However, in the BSG problem the voltage violates the charge 
conservation. Indeed, this is responsible for the charge tunneling in the Luttinger liquid 
with an impurity. In the presence of a voltage in the BSG model, current fiows, so we 
can compute transport properties using a kinetic equation [[7|, p!0| , pT| , |T7| . However, this is 



not useful in finding the free energy at non-zero voltage. Therefore one can use the series 
expansion as above, or as below, we will infer the free energy by utilizing a (well-checked) 
conjecture. As a byproduct, we will also find some more information about the transport 
properties. In particular, we conjecture relations for the conductance good for all values 
of g, generalizing the results of p|, p!0| , pT[1 . 

The partition functions can be expanded in powers of x as before: 



oo 



Zi{x,p) = 2cosp7r -I- y^a:^"Q2n(p) 

n=l 

oo 

ZbSg{x,p) = 1 + Y,X^"l2n{p), 



(4.1) 



71=1 



where 



Qn{p) = J dui...J du'^l2n{{ui},{Ui})2cospin + ^{ui 
^2n(p) = y dui...j du'^X2ni{ui},{u'i})exp iip'^ix 
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where p = igV/2TTT. In |1C] exact series expressions for the l2n{p) were found for integer 
p: 

, , 1 TT r [rrii + g{n - i + l)]r \p + + g{n - i + 1)] 

^2n[P) Y{g)^^ ^i-J^T[m, + l + gin-i)]T\p + m, + l+g{n-i)y ^ ' 

where m is defined as in ( p.6|) . As before, this series converges only for g < 1/2. 

Even though these results apply formally only to p integer, the above formula can be 
applied to all p in the complex plane. To prove that this is the unique analytic continuation, 
we need to assume analyticity as p ^ oo. We know that that /2n(p)/r^"'-^~^^ is analytic 



as p ^ zoo (where T ^ 0) from the analysis of |T^], and we assume this applies at real p 
as well. This assumption has been checked with TBA results below. The series for I2 can 
be summed as before: 



. , >^ T{g + mi)T{g + Xi + p) 
sirnxg r(l - 2g) 



r2(o)r(i + mi)r(i + mi + p) 

mi=0 \ ' ^) \ ' ^ ' 2^ 



sin ^{g + p)V{l - g + p)T{l - g - p) 



One can indeed check that this has the appropriate limit as p ^ zoo to reproduce the 
zero-temperature coefficient of [0, giving support to our analyticity assumption. Notice 
that the continuation of l2n{p) is not even in p, and that Zbsg{p) is I'sal only for p integer. 
The true partition function can be defined by non-equilibrium methods like the Keldysh 
formalism. However, we will see that observable quantities like the conductance are given 
in terms of Zbsg{p)- 

The analysis of sect. 2.2 can be repeated for p non-zero. The fusion relations ( |2.14| ) 
apply without modification for spin-j/2 representations. The /c-dimensional cyclic repre- 
sentation when q'^ = ±1 also obeys the same relation, but to find the BSG free energy, 
there is a subtlety. We first note that when p 7^ we have 



Zs{x,p) ^ I J^e^-^H Zbsg {j^^^'^p) ' q' = C » I (4.4) 



when fundamental set of values is taken to be 0, ... , k—1 (recall that periodic representa- 
tions are invariant under overall shifts oi Sz)- When one fuses such a spin-5 representation 
with a spin 1/2 representation, the last relation in ( |2.14| ) still applies, but the new cyclic 
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representations have shifted fundamental sets, l,...,/c and — — 2, respectively. 
Therefore the relation between Zx and Zbsg is slightly modified: 

Zi[{q-q )x,p] = ^ . (4.5 

Zbsg{x,P) 

Plugging in the power series expansions into gives the hnip) in terms of the Q2n{p)- 
For example, 

Hp) = - — r~^'52(p) 

sm {g + p) 

By direct, straightforward integration, one can check that 

r(l - 2g) 



Q2{P) 



r(l-g + p)T{l-g-p) 



in agreement with (|4.3|) . At next order, we find similarly 

4:sin^{7Tg)Q4{p) = -2 cosTrg sm[n:{2g - p)]l4{p) + sin[7i{g - p)]l2{p)'^ . (4.6) 

generalizing the p = relations in This relation ( [4.5| ) has been checked, again by 

numerical determination of the TBA results for the Kondo model at non-zero magnetic 
field. 

The functional relation ( [4.5| ) implies even more than the Kondo partition function. For 
example, we know on physical grounds (and from the definition ( [4.1D ) that the Kondo par- 
tition function is even in p. Combining this with ( [4.5|) then yields a non-trivial functional 
relation for Zbsg- 

gZTTp lZBSGiqX,p)ZBSGix,-p) - ZBSGiq~^X,-p)ZBSGix,p)] = 

(4.7) 

e '""P [ZBSGiqx, -p)Zbsg{x,p) - ZBSGiq ^x,p)Zbsg{x, -p)] . 

Plugging in the perturbative expansion gives the values of the coefficients l2n{—p) in term 
of the hnip)- 

hip) sm[n{g + p)] = hi-p) sm[n{g - p)] 

h{p) sin[7r(25f -l-p)] - h{-p) sin[7r(25f - p)] = - sin[7rp]/2(p)/2(-p), 
for example. 

In Ijl^, we made and checked a conjecture which related the linear-response conduc- 
tance directly to the partition function. Our conjecture at general fi = V/2T and g is 

^'2T^ ^ 2 d{V/2T)dx \ZBSG{x,-iV/27vtT) J ^ ' 
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We have checked this generahzed expression numerically as well, by comparing it to the 
conductance from the Boltzmann equations in []^,|TT[] and below. 

As at y = 0, much is learned by comparing the perturbative and TBA analyses. The 
TBA equations (|2.16| ) aX g = 1/t are modified slightly in the presence of a finite voltage, 
yielding |,| 

= IZ^jfc *ln(l + e^'=e^'=). (4.9) 
k 

where the incidence diagram is as in sect. 2. A magnetic field corresponds to the chemical 
potentials fi± = ±fi, fik = otherwise and fx = V/2T. The two end nodes of the incidence 
diagram now play a different role, but we still have e_(_ = e_ = et-i- The relation 



still holds. One also has 



liiZt-i{x, fi) = st-i * [In (1 + e'^e^*-!) + In (l + e-^e''-')] 



We then define 



liiZ±{x,n) = --In 



1 + e' 



, smh[{t - l)fi/t] 



+ st_i*ln(l + e^''e'*-i). 



smh[ii/t] 

Using the same identity ( |2.27| ) as for /i = yields 

.s'mh[{t-l)ij/ty~^ 



(4.10) 



(4.11) 



Z±(qV2^,/x)Z±(q-1/2x,/x) 



1 + e 



sinh|U 



sinh[///t] 



(4.12) 



where 



sinh[/i/t] 



Z+Z_. 



Similarly, relation ( |2.31|) becomes 

Zt{x, n) = Zt-2{x,iJt) + 2 cosh//. 



(4.13) 



since the two additional states in the spin t representation have third component of the 
spin equal to ±t. The fusion identites like ( p.29| ) carry over to the case of finite voltage. 
Following the same arguments as for |U = one finds then 



Z_(x,//) 
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(4.14) 



where we remind the reader that fj, = V/2T. 

The functions Z± are not obviously related to any Kondo-type integrals. Since there 
are technical obstacles to directly calculating the boundary sine-Gordon partition function 
in the presence of a voltage, we proceed using the algebraic approach. Comparing ( [4.14| ) 
and (^]^) suggests the functional relations 



together with 



sinh[(U/i] 



Zbsg{x,'^-^)Zbsg{x,-'^-^) . , 

nt Tct smh n 



n-i 



(4.15) 



(4.16) 



Here we traded the p variable of ( [4.1| ) for the fj, variable. It is very likely that ( [4.16| ) 
has an algebraic origin. This is because the tensor product of two cyclic representations 
decomposes on pairs of (generally) indecomposable representations, which in turn are 
related with the spin t — 1 representation of vanishing g-dimension. However, we have 
failed in finding a complete algebraic proof of ( |4.16D , but we have checked it thoroughly by 
using the series expressions above for the Zbsg and the numerical TBA results for Zt-i- 
We can check that ( |4.15D is consistent with the conjectured relation ( [1.8| ) between 
the partition function and the conductance. Using the TBA and a kinetic equation, the 
conductance at integer t = 1/g is p|,pT 



T{t-1) d 



de 



2 dV J cosh\t - l){e - a) 



In 



1 + e-^/2Te- 



■et-l 



(4.17) 



Using the identity 



lim 



it follows that 



1 



cosh (6* + z7r/2 — x) 



cosh [6 — 171/2 + x) 



(4.18) 



-G{q 



' 2T^ 



iiTX d 



d 



2t dxd{V/2T) 



In 



1 + e 



V/2T 



1 + e-v/^'^e'*- 



(4.19) 



This allows a powerful check on the conjectures ( |4.8| ) and ( |4.15| ), because it also follows 
from subsituting ( [4.15| ) into ([4.8|)and using the definition of Z± ( [4.11| ). It would be nice 
to reverse the order of the proof and show that ( [4.19| ) (known to be true from the TBA) 
implies (|4.8| ) and ( |4.15| ). This cannot be done by substituting the perturbative expansion 
because the relation ( [4.19| ) does not determine all of them uniquely; the order x^^ term 
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vanishes on the left-hand side for any integer j. However, it is conceivable that by exploiting 
additional analyticity information that it could be proven along the lines of |28[j . 

In the linear-response limit V ^ 0, we can recover another functional relation from 
T0[] . In this limit we can ignore the V dependence of Yj-i because it is a function of V^. 



Using (|2.30|) , we recover 



Giq'/'x, 0) - Giq-'/'x, 0) = ing'x^- , , (4.20) 

This formula is nice because it no longer has any reference to the TEA quantities e, so we 
expect it to hold for all g. 

We have a formula ( [4.16|) which relates the product of Zbsg{V) and Zbsg{~V) to 
TEA quantities, and a formula ( [4 .81 ) which relates their ratio to TEA quantities by using 
( [4.17I) . Therefore, we can infer a complete expression for Zbsg{x, /u) alone in terms of the 
TEA quantities: 

4 21) 

H nr, . wn ^ In 



where 



e 



1 + g-2(t-l)(0-«) yij^ get_i(-oo)gV^/2T 

et-i(-oo) _ sinh[(t - 



sinh[/x/t] 

and a = ln(TB/T) and fx = V/2T as always; x is related to a via ( p.2Cl|) . It should be 
possible to derive this directly from the TEA, but there are some technical obstacles. 

Extending this analysis to the repulsive regime is more difficult. The reason is that 
for integer s in g = 1 — 1/s, the value of 5"^ at the top and bottom states of the spin-s 
representation still is ±s, not as would be needed to carry over the algebra of the 
attractive regime straighforwardly. This issue is related to the q versus —q~^ problem we 
had to address in sect. 3.2. We find after some manipulation that 

ZBSG{q^^^x,ifi/ns) _ Z+{q^/^x,fi) 



ZBSG{q~^^'^x, ih/tis) Z_ {q-^/'^x, ji) ' 

as above, but with however the new correspondence = (s — 1)^. 
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We can also propose a formula for the conductance in the repulsive regime. By using 
the conjectures (|4.16|) and (|4.8| ) and the identity ( |2.27| ), we now deduce 



ITXX 



l\ d 



d 



sj dx d{V/2T) 



In 



;L _^g(s-l)l//2Tg-e,_i 
1 + f,-{s-l)V/2T(,-es-i 



from which we obtain 



G ( X, 



V_ 
2T 



d9- 



X 



d 



diy/T) 



In 



1 

cosh2(^- In Ts/T) 

1 4- g(s-l)l^/2Tg-e,_i(e) 
1 + e-(s-l)V^/2Tg-e,_i(e) 



In 



1 + g(s-l)V/2Tg-e,_i(oo) 
1 4- e-(s-l)V'/2Tg-£3-i(oo) 



(4.23) 

By generalizing the result of eqn. ( |2.21|) to non-zero V , one can check that these formulas 
give the correct limit G(0, y/2T) = (1 — In the limit of vanishing voltage (linear 

response), one finds 



G{x,Q) 



is -I] 



de 



cosh^ie -IhTb/T) Vl + e^-iW 1 + e 



(oo) 



(4.24) 



This expression has been compared to real-time Monte Carlo simulations in [29|; the 



agreement is good. At T = 0, it agrees with the expression derived in [^. Using the 
identity ( [4.18|) with (|4.24| ) and using the TBA expression for the free energy from sect. 
3.2 yields the functional relation (|4.20 ) in the repulsive regime, lending support to the 
conjecture that ( [4.20| ) holds for all g. 



5. Conclusion 

It should be possible to use the identification of the boundary sine-Gordon model with 
a Kondo type problem more completely than we have done. One way of doing so would be 
to write and solve the Bethe ansatz equations for an integrable system made of spins 1/2 
and a cyclic impurity. Since (|1.3|) conserves the charge, unlike (L2), this would allow one 



to handle directly the BSG model with a voltage, avoiding the lengthy series of functional 
identities of section 4. 

Another interesting direction is to try to continue the perturbative Anderson- Yuval 
coefficients past g = 1 into the irrelevant regime. This can be done by our zeta-function 
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trick of section 3: use the explicit series expression to find out how the series is diverging, 
and subtract and add the appropriate zeta function. This is straightforward but tedious, 
so we have not completed this program. For example, the duality g 1/g should 
be explicitly observable. So far this duality has been established only at vanishing tem- 
perature [0 . Observe however that such a way of handling irrelevant operators does not 
involve any cut-off, and will not describe the non-universal physics depending on the cutoff 
(e.g. the dissipative quantum mechanics in 13^]). However, our result does provide the 
exciting prospect of a well-controlled irrelevant perturbation theory, defined by an analog 
of dimensional regularization. 

Acknowledgements 

We thank S. Chakravarty, R. Egger, M.P.A. Fisher, A. Ludwig, K. Leung, C. Mak, 
J. Rudnick, A. and Al. Zamolodchikov for useful discussions. Some of these results were 
obtained independently by V. Bazhanov, S. Lukyanov and A.B. Zamolodchikov (unpub- 
lished). This work was supported by the Packard Foundation, the National Young Inves- 
tigator program (NSF-PHY-9357207), the DOE (De-FG03-84ER40168) and a Canadian 
NSERC postdoctoral fellowship (F.L.). 



Appendix A. Deriving the S matrix for boundary sine-Gordon 



The boundary S matrix of the boundary sine-Gordon model was found in by 



analyzing the most general solution of the boundary Yang-Baxter equation. Here we show 
that this form also follows from the identification of the BSG model with a cyclic-spin 
anisotropic Kondo model. 

For the ordinary spin-j/2 Kondo model, the S matrix for a particle scattering off the 
impurity is easy to obtain. Up to an overall proportionality factor which follows from 
crossing and unitarity, it is simply the standard Yang-Baxter solution for a spin [j — l)/2 
and a spin 1/2, and a renormalized quantum group parameter q = e*'^^/*^^"^^ [^. By 
analogy, we expect the S matrix for the cyclic spin case to be given (up to the overall 
factor), by the Yang-Baxter solution for a cyclic spin and a spin 1/2. This R matrix is the 
object studied long ago . It is conveniently written as a matrix in Hi C?) H^ 



t^o^'o + 2m;3S'3 2wiS- 

2wiS+ WqSq-2w3S3 



(A.l) 
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where 

a = sin(7 + u),b = sinu, c = sin 7, (A.2) 



and 

a + b a — b 

Wo = , W3 = , wi = 



and in the cychc representation, S± act as indicated above, while 



So\m > 
Ss\m > 



2(gl/2_g-l/2)- 



Setting u = — 6) + v and taking as above the hmit = C » 1, one finds the R 
matrix 

,+,m 



R+,m - e q 

-,m 



(A.3) 



T,+,m-l _ 



m 



The result of [|2^ is that the boundary sine-Gordon S matrix provides a solution of the 



boundary Yang-Baxter equations of the form 

fe-'" 1 

R=[ . I , (A.4) 

while ( [A.3|) is a solution of the ordinary Yang-Baxter equation. To map the two, it is 
tempting to simply forget the cyclic degrees of freedom, which appear only as rapidity- 
independent phases. It is however not totally possible, and this has to do with the difference 
between BYB and YB even for massless particles. Indeed in a massless theory the left- 
right scattering is rapidity independent, but it might still involve some phases. As such, 
( |A.4| ) solves BYB but does not solve YB, because of the left- right scattering phases. 



When considering YB, there are no left-right scattering phases, so the equivalent terms 
are furnished by the cyclic degrees of freedom in (|A.3| ). This is illustrated in figures 4 and 
5. The complete translation shows that, if ( [A.3| ) satisfies YB, then (|A.4| ) satisfies BYB 



indeed. In the figures 4a, 6, c we consider a particular case of the YB equation involving 
the scattering of two spins 1/2 and a "cyclic" spin. The weight of the first figure is Wi = 
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a{u — v)e ^^mt^ ^^le one of the second figure W2 = h{u — v)e *'"e *'^'^/*e ^'^i'm+i)/t^ 
and tfie tliird W3 = ce-^e-^^'^""/*. Tlie fact that YB holds means that Wi = W2 + W3. 

In the figures 5a, 6, c we consider a particular case of the BYB equation involving a 
spin 1/2 bouncing off a boundary. The weight of the first figure is W{ = e~^'^aLL{u — 
v)aLR{u — f), the weight of the second figure W2 = e~^'^bLL{u — v)bLR{u — f), and the 
weight of the third = e~^'"cRR{u — v)aLR{u — v). The fact that BYB holds means 
W[ = W2 + VF3, which one checks easily since the RR elements are identical to the ones 
in and ulr = 1, bLR =^ e"'^, clr = 0. 
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Fig. 4a: Configurations involved in checking YB equation with two spin 
1 /2 and a cychc representation (see the appendix) 
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Fig. 5a: Configurations involved in checking BYB with no degree of freedom 
at the boundary (see the apppendix) 
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Fig. 5b, c: idem 



40 



